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Newton’s Law of Cooling: An Introduction to Scaling 
1 Introduction 


This simulation applies scaling to solve a relatively simple differential 
equation involving Newton’s law of cooling. Initially at 100°C, a hot cup of tea 
is placed in a room that is maintained at an ambient constant temperature of 
20°C. 


For comparison, analytical, numerical, and analog computer solutions are 
displayed in table 1 on page 6. 


Please note: Originally, the rocket equation derivation was going to be part of 


Application Note #9. Instead. it will be included in Application Note #10. 
Sorry for any inconvenience this may have caused. 


2 Mathematical modeling 


Starting with Newton’s law of cooling, 
dT/dt = -k(T - Ta) with T(0) = To, where (1) 
T = temperature of the tea 
t= time 
k = temperature coefficient 


Ta = ambient temperature (assumed constant) 
To = initial temperature of the tea 


Letting k = 0.25 min-!, Ta = 20°C, and T(0) = 100°C, 
dT/dt = -0.25min“+(T — 20°C) with T(0) = 100°C (2) 


Following a bit of calculus (details will be provided for the scaled version), 
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T(t) = 80°Ce*(-0.25 min't) + 20°C (3) 
As a reminder, operational amplifier (op amp) input/output voltages are 
between Vez and Vcc. For this project, 9-Volt batteries were used (Vez = -9 
Volts and Vcc = +9 Volts). Allowing for a safety of margin, to prevent 
saturation, voltages are kept between -6 Volts and +6 Volts. 
The issue of magnitudes is obvious. A value like 100°C is just too big for direct 
conversion to 100 Volts. So, a temperature of 100°C will be reduced (scaled 
down) to an analog voltage of 6 Volts and 20°C will be reduced to an analog 
voltage of 1.20 Volts. 
Time will be scaled such that 1 minute = 1 second. Easy enough! 
To start, let 
R = Reducing scale factor = Tmax/Vmax = 100°C/6 Volts = 50°C/3 Volts. 
In general, R= T/V or T= RV. 
Replacing T with RV, (1) becomes 
d(RV)/dt = -k(RV - Ta) with V(0) = T(0)/R = To/R 
RdV/dt = -kR(V - Ta/R) with V(0) = T(0)/R= T./R 
dV/dt = -k(V - Va) with Va= Ta/R and V(0) = T(0)/R=T./R 
Inserting values, but omitting units for clarity, 
dV/dt = -0.25(V - 20/(50/3)), V(0) = 100/(50/3) 
dV/dt = -0.25(V - 1.20) with V(0) = 6.00 (4) 
dV/(V - 1.20) = -0.25 dt 
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Integrating by inspection, and noting that V > 1.20, 
In((V - 1.20)/(6 - 1.20)) = -0.25t 
In((V - 1.20) /4.8) = -0.25t 
(V - 1.20) /4.8 = e*(-0.25t) 
V = 4.80 e*(-0.25t) + 1.20 
Inserting units, 
V(t) = 4.80 Volts e*(-0.25 mint) + 1.20 Volts (5) 
Since T = RV, 
T(t) = 50°/3 Volts x [4.80 Volts e*(-0.25 min“!t) + 1.20 Volts] 
T(t) = 80°C e*(-0.25 mint) + 20°C, which is identical to (3) 
3a Computer setup (scaled, patch cord version) 
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Figure 1: Computer setup for Newton’s law of cooling 


3b Computer setup (op amp/discrete component version) 
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Figure 2: Basic breadboard layout 


4 Numerical Method 
(Modified Euler method using a hand-held programmable calculator) 
Code: TI-BASIC 


PROGRAM:COOLING 
:ClrHome:ClrDraw 
:’NEWTONS LAW” 
"OF COOLING” 
:"DV/DT=-0.25(V—1.20)” 
:"WITH V(0)=6.00” 
:’V=VOLTAGE ANALOG” 
?”?OF TEMPERATURE” 
?”T =TIME” 
:”’PARAMETERS:” 
:0-T:6.00-V:0.25-H 
:”"H=STEP SIZE” 

: Fix 1 

:Lbl 1 

‘If T>20:Then 

:Goto 2: Else 

:DisP {T,V} 
:-0.25(V—1.20)>F 
:V+HF>W 

:‘T+HOT 
:-0.25(W-1.20)-S 
((F+S)/2-A 
:V+AH-V 

:Pause 

:Gotol 

:Lbl 2 

‘End 


5 Results 
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Figure 3: Voltage vs time simulation* 


*For this application note, the oscilloscope display was produced during a 
single run using a differential equation analog computer constructed from 
operational amplifiers and discrete components with tolerances within 10%. 



































t (s) scaled Analog Analog Analytical Numerical 
from minutes Computer Computer T (°C) T (°C) 

V (Volts) Converted 

Estimated from | T (°C) 

oscilloscope 
00.0 6.0 100 100 100 
02.5 3.6 60 63 63 
05.0 2.6 43 43 43 
07.5 2.0 33 32 32 
10.0 1.6 2/ 27 zy 
12.5 1.3 ZZ 24 24 
15.0 12 20 22 22 
17.5 1.2 20 Z1 ZL 
20.0 1.2 20 21 21 




















Table 1: Solution Comparisons 


6LP3146 MNI604 MAR2024 
6LP3146 MNI604 MAR2O23 


DURACELL? 


9V 


DURACELL’ 


9V 


a! 1.5A 


WERED BREADBOARD 
MODEL:PBB-4060B 


(oe wee 
eee eee 














Figure 4: Differential Equation Analog Computer 
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